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Abstract 

Temperature  and  hydration  management  are  very  important  points  for  operating  fuel  cell  systems  and  for  optimizing  their  energy  balance. 
In  this  paper,  the  design  and  the  manufacture  of  an  air  conditioning  device  dedicated  to  proton  exchange  membrane  (PEM)  fuel  cells  are 
presented.  Experimental  thermal  behavior  laws  are  described  and  are  intended  to  be  integrated  in  a  fuel  cell  powertrain  model  to  drive  the  fuel 
cell  without  using  an  expensive  and  not  adapted  hygrometry  measurement  system.  Experiments  on  a  low  power-testing  bench  (500  W)  have 
been  carried  out,  showing  the  efficiency  of  the  system.  Fuel  cell  polarization  curves  are  shown,  overheating  and  thermal  inertia  durations  are 
also  provided  and  taken  into  account  for  optimizing  the  regulation. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Hydration  content  of  the  fuel  cell  membrane  is  a  key  point 
to  get  a  good  efficiency  [1,2]  of  a  PEM  fuel  cell  stack.  The 
produced  water  during  electrochemical  reaction  at  the  cath¬ 
ode  side  is  not  usually  sufficient  to  guarantee  good  moisture 
content,  in  the  whole  range  of  operating  current  density  and 
temperature.  Additional  water  has  to  be  brought  to  the  system 
by  hydrogen,  air  or  both  [3].  The  water  balance  is  not  easy 
to  reach  as  too  much  water  leads  to  a  flood  of  the  electrodes 
and  not  enough  water  will  dry  the  membrane  [4] .  The  behav¬ 
ior  of  the  stack  versus  this  phenomenon  depends  on  a  great 
amount  of  parameters:  type  and  thickness  of  the  membrane, 
catalyst  layer,  electrode  backing  and  flow  channels,  temper¬ 
atures  of  the  different  parts  of  the  stack.  Furthermore,  the 
amount  of  water  dragged  by  the  gas  depends  on  the  humidity 
rate  of  the  gas  and  on  its  temperature.  Moreover,  mastering 
the  gas  temperature  and  its  humidity  rate  costs  a  lot  of  en¬ 
ergy  and  lowers  the  global  efficiency  of  the  fuel  cell  system. 
In  order  to  study  this  critical  issue,  an  air  conditioning  sys- 
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tern  has  been  designed  and  realized  for  a  low  power  PEM 
fuel  cell  test  bench  (500  W).  After  having  described  different 
hydration  systems  and  justified  our  choice,  the  second  part 
of  this  paper  describes  the  testing  bench  [5,6]  and  the  prin¬ 
ciple  of  the  air  conditioning  system.  The  choice  of  the  boiler 
for  this  application  is  valid  for  obtaining  a  reliable,  dynamic 
and  low  cost  hydration  management  system,  moreover  easy 
to  control.  It  has  been  built  and  adapted  to  the  low  power 
of  the  PEMFC.  The  following  part  of  the  paper  shows  ex¬ 
perimental  results  to  demonstrate  the  ability  of  the  system  to 
ensure  both  a  given  humidity  rate  and  a  given  temperature  of 
the  input  air.  The  thermal  inertia  of  the  boiler  and  the  heater 
are  also  shown.  The  last  part  describes  thermal  behavior  ex¬ 
perimental  laws  of  the  air  conditioning  system.  The  thermal 
characterization  laws  are  obtained,  thanks  to  experimental 
temperatures  and  time  durations  of  the  air  conditioning  an- 
cillaries  transient  states.  Theoretical  laws  permitting  to  have 
enthalpies  as  function  of  air  pressure  are  also  calculated.  Law 
association  is  then  integrated  in  an  open  loop  control  of  the 
air  conditioning  device.  It  finally  permits  to  control  the  air 
temperature  and  hygrometry  levels  without  using  a  very  pre¬ 
cise  (and  also  very  expensive)  optical  hygrometer  that  cannot 
be  considered  for  PEM  fuel  cell  embedded  applications. 
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Nomenclature 

A,  B,  C,D  constant  coefficient  (no  dimension) 

CP 

specific  heat  (Jkg”1  K”1) 

Cp 

heat  capacity  (J  K”1) 

h 

specific  enthalpy  (Jkg”1) 

h'm 

input  air  specific  enthalpy  (J  kg”1 ) 

^out 

output  air  specific  enthalpy  (J  kg” 1 ) 

hr 

saturated  air  specific  enthalpy  (J  kg” 1 ) 

Pi 

absolute  pressure  (Pa) 

Pv  s 

saturated  air  absolute  pressure  (Pa) 

P  ext 

external  heat  losses  (W) 

Pv 

feeding  power  (steady  state)  (W) 

Pf 

feeding  power  (transient  state)  (W) 

Pi  hi -2 

thermal  power  transmitted  to  the  air  (W) 

dm 

mass  flow  (gl”1) 

qv 

volume  flow  rate  (N 1  min” 1 ) 

Q 

amount  of  heat  (J) 

T 

dry  air  temperature  (°C) 

Text 

room  temperature  (°C) 

Tair  input? 

7jn  input  temperature  of  the  air  (°C) 

Tair  outpu 

t,  Tout  output  temperature  of  the  air  (°C) 

Tr 

saturated  air  temperature  (°C) 

Greek  symbols 

At 

time  duration  (s) 

P 

density  (kgm”3) 

X 

moisture  content  (kgHoQ  kg”1) 

2.  Different  kinds  of  hydration  management  systems 

This  section  is  devoted  to  a  short  state  of  existing  hydration 
systems  already  available  for  PEM  fuel  cells. 


2.1.  Water  injection 

Heated  water  is  injected  in  a  column  where  air  circulates. 
A  controlled  airflow  enters  at  the  bottom  of  the  column  while 
a  controlled  water  stream  is  injected  at  the  top.  As  the  air 
rises,  it  gets  wet  [7].  Vapour  exchange  is  favoured  by  metal¬ 
lic  rolls  inside  the  column.  The  water  quantity  to  be  injected 
is  computed  versus  the  air  temperature,  pressure,  flow  and 
the  required  humidity  rate.  A  bleeder  located  at  the  column 
bottom  collects  water,  which  has  not  been  transferred  to  the 
air.  Water  is  heated  before  entering  the  system  by  an  evapo¬ 
rator  and  its  temperature  is  one  of  the  parameters  to  be  con¬ 
trolled.  However,  the  exchange  efficiency  is  less  than  1  as 
liquid  water  in  the  bleeder  witnesses.  As  this  amount  of  liq¬ 
uid  water  is  difficult  to  evaluate  accurately,  there  is  a  discrep¬ 
ancy  between  the  actual  humidity  rate  and  the  humidity  set 
point. 


Warm  Humid 


Inlet  Air 

Fig.  1.  Enthalpy  wheel  humidifier  [8]. 


2.2.  Enthalpy  wheel  humidifier 

This  device  shown  in  Fig.  1  reuses  the  enthalpy  of  the  out¬ 
let  gas  to  humidify  and  warm  the  inlet  fresh  air.  It  consists 
in  a  porous  cylindrical  wheel,  coated  with  desiccant,  rotat¬ 
ing  slowly  inside  the  housing.  Warm,  humid  exhaust  air  from 
the  fuel  cell  enters  the  humidifier.  The  rotating  drum  absorbs 
the  exhaust  water  vapour  transferring  it  to  the  cool  dry  in¬ 
let  stream.  Other  exhaust  gases  continue,  thanks  to  the  drum 
desiccant  coating.  This  device  simplifies  greatly  the  system, 
it  avoids  the  output  condenser,  and  it  saves  energy  as  it  recov¬ 
ers  exhaust  enthalpy.  However,  it  is  hard  to  control  accurately 
the  humidity  and  the  temperature  of  the  inlet  air.  It  depends 
on  the  temperature  and  humidity  of  the  exhaust  air,  i.e.  on  the 
stack  operating  point.  It  might  also  introduce  instabilities  in 
the  system  control.  This  attractive  solution  can  be  useful  in 
an  operating  system,  as  long  as  the  fuel  cell  stack  has  been 
deeply  characterized  versus  the  air  parameters  and  the  hu¬ 
midifier  performances  are  well  known.  It  should  be  firstly 
studied  and  compared  with  a  laboratory  device,  which  al¬ 
lows  repeatable  air  conditioning,  independent  from  the  stack 
loading. 

2.3.  Boiler 

This  device  consists  in  the  injection  of  an  external  air  in 
the  bottom  of  a  bottle  filled  with  water.  It  is  the  system  that  is 
used  in  the  first  stage  of  the  proposed  air  climate  control  unit. 
The  required  temperature  of  the  water  and  finally  of  the  air 
obtained  at  the  output  of  the  boiler  is  ensured,  thanks  to  ther¬ 
mal  resistances.  The  required  output  air  hygrometry  of  such 
ancillaries  leads  to  100%,  whereas  a  following  heater  allows 
obtaining  different  values  of  air  temperature  and  hygrometry 
at  the  entrance  of  the  stack.  The  initial  choice  of  the  device  is 
accounted  for  obtaining  a  reliable,  faithful  and  not  expensive 
hydration  system  allowing  minimizing  the  thermal  transient 
duration.  This  point  is  a  very  important  one  considering  fu¬ 
ture  dynamic  applications.  It  also  permits  to  work  in  a  very 
large  temperature  and  hygrometry  band  and  to  have  a  stable 
regulation  loop. 
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3.  Testing  bench  description 

3.1.  The  stack 

The  stack  under  test  has  been  supplied  by  ZSW  Germany. 
It  is  composed  of  20  cells.  The  active  area  is  100  cm2.  It  is  fed 
on  the  anode  side  with  dry,  pure  hydrogen.  It  can  be  fed  on 
the  cathode  side  either  with  pure  oxygen  or  air,  dry  or  humid. 
The  operating  temperatures  take  range  from  20  to  70  °C.  The 
relative  gas  pressure  has  to  be  kept  under  0.5  bar  and  the 
pressure  drop  between  anode  and  cathode  must  remain  under 
250  mbar.  The  cooling  water,  circulating  in  the  bipolar  plates 
must  be  deionized  (electrical  conductivity  <2  jxS  cm-1).  The 
maximal  electrical  power  output  is  600  W  under  12  V  [9]. 

3.2.  The  air  conditioning  design 

This  paper  highlights  the  design  of  the  air  circuit  for 
the  control  of  both  its  temperature  and  humidity  rate.  The 
operating  range  of  the  air  conditioning  is  bounded  by 
two  equal  humidity  rate  curves:  xi  =0.015  kgvap  kg and 

X2  =  o.  1 5  kgvap  kggaj, . 

Fig.  2  shows  the  air  conditioning  circuit  principle.  It  is 
made  of  four  elements: 


-  the  boiler  saturates  air  at  the  chosen  due  point  temperature; 

-  the  condenser  cools  the  saturated  air  to  reach  the  desired 
amount  of  water; 

-  the  separator  evacuates  the  condensed  water; 

-  the  heater  warms  the  gas  to  the  required  temperature  and 
due  point. 

The  airflow  is  determined  according  to  the  load  current 
and  controlled  through  a  mass  flow  regulator.  The  due  tem¬ 
perature  is  calculated  according  to  the  required  temperature 
and  hygrometry,  and  the  corresponding  absolute  pressure. 
Dry  (or  not)  incoming  air  goes  firstly  through  the  boiler.  At 
the  boiler  output,  air  presents  a  100%  hygrometry  level  and  a 
temperature  Tr.  A  discrepancy  between  the  due  temperature 
Tv  and  the  dry  temperature  set  point  T  implies  the  warming 
of  the  boiler.  Then,  wet  air  is  cooled  in  the  condenser  to  ex¬ 
tract  the  required  amount  of  water  from  the  air  (in  order  to 
obtain  the  required  air  hygrometry  level  at  the  stack  input). 
The  separator  ensures  the  evacuation  of  the  condensed  water. 

Finally,  the  temperature  is  also  measured  at  the  stack  input 
(Fig.  2).  In  the  case  of  a  difference  between  the  measured 
temperature  and  the  required  input  gas  temperature,  the  heater 
can  be  switched  on. 

In  order  to  check  the  reliability  and  faithfulness  of  the 
air  management  device,  an  optical  hygrometer  has  been 


Fig.  2.  Air  conditioning  system. 


Dew  point  (Tsat)  and  inlet  Air  temperature  (TeAir)  behaviour 


Fig.  3.  Testing  bench  temperature  responses. 
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Polarization  curves  with  various  input  air  dew  points  Tair 
stack  temperature  T  =  50  °C 


Fig.  4.  Polarization  curves  measured  under  different  dew  points  of  the  input  air. 


temporarily  placed  at  the  stack  input.  However,  it  is  not  fi¬ 
nally  involved  in  the  process  regulation  because  of  its  low 
time  dynamic  response  (shown  in  Fig.  3)  and  its  high  price. 
The  whole  hygrometry  setting  system  is  also  operating  in  a 
so-called  open  loop  mode. 

4.  Experimental  results 

The  air  conditioning  device  allows  investigating  the  in¬ 
fluence  of  the  temperature  of  the  input  air  and  its  humidity 
on  the  stack  behaviour.  The  aim  of  this  type  of  study  is  to 
determine  optimized  operating  conditions  for  the  whole  fuel 
cell  system.  The  obtained  temperatures  and  hygrometry  in¬ 
ertia  and  reliability  are  also  shown  in  this  sectionmpared  to 
the  required  ones. 

4.1.  Temperatures  reliability  for  a  100%  required 
humidity  rate 

A  test  has  been  realized  for  determining  temperature  reli¬ 
ability  for  a  100%  humidity  rate  at  the  stack  air  input.  Three 
operating  points  at  the  dew  point  are  required: 

Tv  =  25  °C,  Fair_  ^pm  =  26  °C/  Tv  =  35  °C, 

Fair-input  =  36  °C/  Tr  =  45  °C,  Fair_input  =  46  °C 

Fuel  cell  power  generation  is  490  W  and  stoichiometry 
factors  at  the  anode  and  the  cathode  are,  respectively,  2 
and  4.  Hydrogen  volume  flows  at  the  input  and  output  of 
the  anode  are,  respectively,  9.86  and  4.78  N1  min-1.  At  the 
cathode,  air  volume  flows  are  in  the  same  order:  46.86  and 
43.43  N1  min-1. 

Fig.  3  gives  temperature  responses  as  function  of  time  t. 
A  relative  reliability  between  requested  and  obtained  tem¬ 
peratures  can  be  seen.  Two  degreescelsius  temperature  dif¬ 
ferences  can  nevertheless  be  seen  between  output  required 
temperature  and  output  obtained  temperature  at  the  boiler. 
This  difference  is  probably  bounded  by  thermal  heat  losses 
of  the  boiler.  Moreover,  the  final  hygrometry  rate  measured  by 
the  optical  hygrometer  is  97%  showing  the  well  behaviour  of 
the  proposed  air  conditioning  system.  Problem  of  Tsat  given 
by  the  hygrometer  is  the  high  thermal  inertia  of  this  device 


compared  to  required  parameters.  Moreover,  such  apparatus 
depends  on  local  thermal  and  fluidic  phenomenon.  Inertia 
and  fall  of  indicate  temperature  between  1200  and  1700  s  are 
bounded  by  local  condensation  drop. 

Generally,  temperature  differences  between  required  and 
obtained  temperatures  come  from  the  thermal  heat  losses, 
particularly  at  the  boiler  and  the  heater.  To  minimize  these 
differences,  it  is  also  necessary  to  quantify  very  precisely 
all  these  losses  at  different  operating  levels.  A  good  char¬ 
acterization  of  these  points  also  permits  to  regulate  the  air 
management  device  faithfully. 

Fig.  4  shows  experimental  polarization  curves  for  a  100% 
humidity  rate.  They  represent  the  stack  voltage  as  functions 
of  the  current  for  different  due  points  temperatures  for  the 
input  air  taking  range  from  22.2  to  43.5  °C.  The  stack  tem¬ 
perature  is  kept  constant  (at  50  °C).  The  stack  output  cooling 
water  temperature  has  been  considered  as  representative  of 
this  stack  internal  temperature.  The  air  enters  the  stack  at  its 
dew  point,  which  varies  between  the  different  curves.  The  in¬ 
fluence  of  this  parameter  on  the  stack  response  is  here  clearly 
demonstrated,  as  the  available  electrical  power  is  increasing 
with  the  inlet  air  dew  point. 

4.2.  Polarization  characteristics — validation  of  the 
device 

Fig.  5  shows  a  polarization  curve,  supplying  the  stack  with 
unsaturated  air  (Fstack  =  50  °C,  Tair  input  =  30  °C  and  hygrom¬ 
etry  rates  are  50%  and  75%).  Software  developed  with  Lab- 
view  permits  to  drive  the  boiler,  heat  exchanger  and  final 
heater  in  order  to  obtain  the  requested  input  fluid  parame¬ 
ters  for  the  stack.  The  software  is  validated  thanks  to  the 
responses  of  the  temperature  and  hygrometer  sensors  located 
at  the  input  of  the  stack. 

5.  Thermal  characterization  of  the  experimental  air 
circuit 

The  proposed  experimental  characterizations  also  permit 
to  define  with  a  very  good  accuracy  the  thermophysical  char¬ 
acteristics  of  the  thermal  losses  of  the  boiler  and  the  heater. 
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Polarization  curve  (Tstack=50°C) 


Hygrometry  level  =  50% 

11  4 
10 


<0  'b  <0  <b  Jo  Q>  Jo  <1  Jo  <0  r<0  <b  r*3  oN  d?  O'  ef°  rS 

N-  V  A-  N  <b*  N  N  <1  fy  ^  ^  cfb-  <b 


qy  *  qy 


Stack  current  (A) 

Fig.  5.  Stack  polarization  curves  with  air  at  50%  and  75%  humidity  rates  (lair  input  =  30  °C). 


Such  identification  also  permits  to  define  precisely  the  be¬ 
haviour  laws  of  the  air  conditioning  system.  Theoretical  laws 
of  the  air  behaviour  [10]  have  permitted  to  obtain  real  ther¬ 
mal  powers  to  bring  to  the  boiler  and  to  the  heater  to  have  the 
required  air  humidity  and  temperature  at  the  input  of  the  fuel 
cell  for  different  fluid  molar  flows  and  absolute  pressures  tak¬ 
ing  range  from  1  to  4  bar.  Then,  experimental  thermal  losses 
of  the  boiler  and  the  heater  are  expressed  as  functions  of  the 
input  and  output  temperatures  of  both  apparatus  for  different 
fluid  molar  flows.  The  sum  of  theoretical  power  and  exper¬ 
imental  thermal  losses  of  apparatus  also  allows  evaluating 
electrical  power  sources  at  the  boiler  and  the  heater  as  func¬ 
tion  of  the  required  humidity  and  temperature  at  the  input  of 
the  fuel  cell. 

5.7.  Air  enthalpy  as  a  function  of  temperature  and 
absolute  pressure 

The  feeding  of  the  stack  is  made  with  atmospheric  pressure 
air.  The  enthalpy  hv  of  the  saturated  air  has  to  be  calculated  as 
a  function  of  the  temperature  T  and  of  the  absolute  pressure  p t. 
An  approximated  expression  of  the  pressure  of  the  saturated 
air  is  defined  by  the  parabolic  equation  (1),  valid  between  10 
and  70  °C. 

pvs  =  9.35T?  -  326Tr  +  5432  (1) 

Enthalpy  hv  is  given  by  equation  (2)  where  x  is  the  humidity 
ratio  and  Tv  is  the  saturation  temperature.  This  expression  is 
given  for  1  kg  of  dry  air. 

hr  =  Tr-\-  x(2490  +  1 .9677)  (2) 

The  humidity  ratio  is  calculated  according  to: 

0.622j?vs 

X  =  -  (3) 


Equations  (1)  and  (3)  in  equation  (2)  allow  obtaining  a 
new  expression  of  enthalpy  given  by  equation  (4): 

_  2.057 T)  +  14402 T}  +  ( pt  -  503452 )TV  +  8408736 
r  “  pi  -  (9.35rr2  -  326Tr  +  5432) 

(4) 

The  evolutions  of  enthalpy  in  a  50-600  kJ  kg-1  range  are 
the  considered  functions  of  the  saturation  temperature  for 
absolute  pressures  from  105  to  4  x  105  Pa. 

A  second-order  polynomial  equation  of  the  enthalpy  as 
function  of  the  saturation  temperature  at  the  considered  ab¬ 
solute  pressures  can  be  easily  obtained  (5). 

h,  =  ATt2  -  BTX  +  C  (5) 

Values  of  constants  A-C  are  given  in  Table  1 .  They  also 
depend  on  absolute  air  temperature  in  the  air  channel. 

5.2.  Thermal  characterization  of  the  boiler  and  heater 
without  external  heat  losses 

5.2.1.  Power  supply  at  the  boiler 

Case  of  a  thermal  steady  state  and  an  absolute  air  pressure 
of  105  Pa,  the  theoretical  power  of  the  system  without  heat 
losses  is  expressed  through  the  balance  of  the  enthalpies  on 
the  entrance  and  output  fluid  channels  (equation  (6)). 

Ahl  —  out  —  ^in)  (6) 

Table  1 


Polynomial  constants  A-C  for  different  absolute  air  pressures 


Absolute 
pressure  /?t  (Pa) 

A  (kJkg-1  °C“2) 

B(kJ  kg"1  “C”1) 

C(kJ  kg-1) 

in 

O 

H 

X 

H 

0.230 

8.43 

133.3 

2  x  105 

0.148 

8.09 

171.1 

3  x  105 

0.085 

3.78 

92.24 

4  x  105 

0.067 

3.34 

92.42 

Pt  Pv  s 
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In  that  expression,  enthalpies  are  calculated  with  equation 
(4).  The  theoretical  power  of  the  boiler  can  be  expressed  as  a 
function  of  the  entrance  fluid  hygrometry  rate  (p  for  different 
values  of  the  required  air  temperatures  T  (25-55  °C).  The 
absolute  pressure  pt  is  105  Pa,  whereas  the  heat  flow  is  about 
qv  =  15  N 1  min-1  and  Tm  =  20  °C.  Software  then  allows  to 
express  the  power  supply  Pthi  as  a  function  of  the  hygrometry 
rate  at  the  atmospheric  pressure  pt.  It  can  be  expressed  thanks 
to  equation  (7): 

Pthi  =  A(T)<p  +  B(T)  (7) 

Expressions  of  A(T)  and  B(T)  can  be  easily  literally  de¬ 
termined  after  having  drawn  their  evolution.  Then,  approx¬ 
imation  laws  of  parameters  A  and  B  expressed  in  equation 
(7)  give  the  final  expression  (8)  of  the  supply  power  Pthi  as 
function  of  required  cp  and  T: 

Pthi  =  (0.0744P2  -  3.4218 T  +  54.397)<p  -  0.0006T3 

+  0.061 1 T2  -  1.74497+  16.373  (8) 

For  a  complete  characterization  of  the  system  in  the  op¬ 
erating  range,  it  is  necessary  to  consider  different  airflow 
alimentations.  Indeed,  different  stoichiometry  factors  can  be 
required  imposing  also  variable  airflows.  Tests  are  realized 
for  different  mass  airflows  available  from  a  minimum  value  to 
1 .6  g  s_1 .  Obtained  approximation  laws  are  given  in  equation 

(9): 

—  =  (0.2978T2  -  13.701  T  +  217.82)^  -  0.0458T2 

+  4.1286T- 69.416  (9) 

All  these  laws  have  been  defined  considering  operating  air 
conditions  at  7in  =  20°C  for  different  airflows  feeding.  The 
temperature  is  the  room  temperature.  Nevertheless,  future 
tests  taking  into  account  variable  external  conditions  will  be 
carried  out  for  higher  PEMFC  powers  thanks  to  a  climate 
chamber. 

5.2.2.  Power  supply  to  the  heater 

In  the  same  way  as  the  boiler,  the  expression  of  the  supply 
power  for  the  heater  is  given  by  expression  (10): 

Pth2  =  (0.0225 T2  -  1.2604 T+  14.204)<p2 

+  (— 0.0156T2  +  0.3033 T  +  4.947)<p  -  0.0025T2 
+  0.6242P- 13.128  (10) 

and  for  different  mass  air  flows,  it  becomes  for  input  condi¬ 
tions  at  pi  =  105  Pa  (1 1): 

—  =  (0.090 IT2  -  5.0413T  +  56.813)<p2 

+  (—0.0624 T2  +  1.2133 T  -  19.789)^ 

-0.0064r2  +  2.1943T- 46.987  (11) 


In  expression  (11),  cp  and  T  are,  respectively,  the  required 
hygrometry  and  temperature  at  the  air  input  of  the  stack. 

5.3.  Thermal  characterization  of  the  boiler  and  heater 
external  heat  losses 

Considering  the  exhaustive  heat  power  generation  at  the 
boiler  and  the  heater,  taking  into  account  conduction,  con¬ 
vection  and  radiation  heat  losses,  only  an  experimental  study 
allows  to  obtain  the  real  thermal  losses  for  the  steady-state 
thermal  heat  transfer.  Obtained  values  of  the  external  heat 
losses  are  to  be  added  to  theoretical  power  supplies  for  each 
element  to  calculate  final  efficiency  of  the  complete  system. 
It  can  be  expressed  thanks  to  equation  (12): 

Text  =  P?  —  Pth  (12) 

with  PY  the  power  actually  supplied  to  the  element,  to  ensure 
output  thermodynamic  constraints. 

Experimental  tests  permit  to  obtain  an  approximation  law 
of  Pext  for  different  required  hygrometries,  output  air  temper¬ 
atures  and  mass  airflows.  The  obtained  laws  are  linear  and 
given  by  equations  (13)  and  (14): 

Textl  —  2. 32(Tair  output  —  Text)  (13) 

Pext2  =  0.1725(Tair  output  Tex  ,)  (14) 

5.4.  Determination  of  the  specific  energy  per  degree  of 
the  boiler  and  the  heater 

Dynamic  tests  showing  thermal  transient  states  of  these 
elements  allow  obtaining  the  specific  energy  per  degree.  Its 
value  translates  the  duration  of  the  transient  thermal  trans¬ 
fer  between  power  generated  and  loose  in  the  boiler  and  the 
heater.  Finally,  it  gives  an  indication  about  the  dynamic  ther¬ 
mal  response  of  the  air  circuit. 

The  expression  of  the  energy  variation  of  the  considered 
system  is  estimated  with  a  transient  experimental  test,  sup¬ 
plying  the  element  with  the  power  PYr  during  a  temperature 
increasing  time  At.  It  is  given  thanks  to  equation  (15): 

Cp(T0Ut  —  Text)  —  [Tr  —  (p^v(^airm  —  fon)  T  Text)]  Af 

(15) 

Parameter  ha irm  represents  the  specific  enthalpy  of 
the  air  calculated  to  the  arithmetic  average  temperature 
(Tout  +  Tjn)/2.  It  is  then  possible  to  obtain  the  specific  en¬ 
ergy  per  degree  Cp  of  the  boiler.  Fig.  6  shows  the  evolu¬ 
tion  of  the  air  output  temperature  during  the  transient  ther¬ 
mal  phase.  Characteristics  of  that  test  are:  Tairouput  =  61  °C, 
Taa  input  =  24.2  °C,  hrm=l  89.24  kJ  kg”1 ,  hm  =  24.2  J  kg” 1 , 
qv  =  0.017  Is.  The  estimated  duration  of  the  thermal  tran- 
sient  phase  is  At  =  1104s.  Calculated  Cp  for  the  boiler  is 
about  28,600  JK-1.  Different  tests  have  shown  a  thermal 
transient  duration  constant  with  t  in  our  operating  range. 

As  shown  in  Fig.  7,  test  characteristics  for  the  determina¬ 
tion  of  the  specific  heat  of  the  heater  are  the  following  ones: 
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Fig.  6.  Temperature  evolution  of  the  boiler  (Ta ir  output  =  60  °C). 


Time  (s) 


Fig.  7.  Temperature  evolution  of  the  heater  (7^  output  =  60  °C). 


^air  output  —  65  °C,  Text  =  22  °C,  heat  generation  Q-  12  kJ  and 
the  initial  temperature  Tm  =  18.6  °C.  The  estimated  duration 
of  the  thermal  transient  phase  is  At  =  1 104  s.  Calculated  Cp 
for  the  heater  is  about  236  J  K_1 . 


6.  Conclusion 

Water  management  in  a  PEM  fuel  cell  stack  is  a  critical 
point.  Therefore,  the  conditioning  of  the  input  air  must  be 
controlled.  The  designed  device,  which  has  been  presented 
and  realized  in  this  study,  allows  controlling  both  the  temper¬ 
ature  and  the  humidity  rate  of  a  PEM  fuel  cell  incoming  air. 
The  influence  of  these  parameters  on  the  stack  performances 
can  thus  be  easily  investigated  on  a  PEM  testing  bench  [11]. 


The  results  allow  defining  the  best  operating  conditions  for 
the  stack.  Moreover,  a  balance  of  the  auxiliaries  consump¬ 
tions  will  lead  to  the  evaluation  of  the  output  net  power  of 
the  fuel  cell  generator. 

A  model  of  the  subsystems’  behaviour  laws  has  been  elab¬ 
orated.  It  allows  determining  the  temperature  references  to 
reach  the  required  characteristics  for  the  fuel  cell  incoming 
air,  taking  into  account  thermal  losses  of  the  boiler  and  the 
heater.  It  also  allows  using  the  hydration  system  without  a 
very  expensive  optical  hygrometer  having  a  too  important 
measurement  time  response.  It  is  then  necessary  to  define 
functioning  laws  of  the  air  circuit,  the  boiler  and  the  heater 
for  different  mass  flows,  air  temperatures  and  hygrometry 
rates.  This  model  is  also  integrated  in  the  open  loop  control 
of  the  air  feeding.  The  described  methodology  applied  to  a 
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specific  kind  of  fuel  cell  is  a  general  methodology  to  simplify 
in  a  near  future  higher  power  fuel  cell  hydration  systems. 
Only  the  constants  of  equations  in  this  paper  will  have  to  be 
adapted  to  higher  power  stacks. 
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